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Competition with magnetism is at the heart of high temperature superconductivity, most intensely felt near 
a vortex core. To investigate vortex magnetism we have developed a spatially resolved probe using nuclear 
magnetic resonance. Our spin-lattice-relaxation spectroscopy is spatially resolved both within a conduction 
plane as well as from one plane to another. With this approach we have found a spin-density wave associated 
with the vortex core in Bi2Sr2CaCu20g+« , which is expected from scanning tunneling microscope observations 
of "checkerboard" patterns in the local density of electronic states. 1 1 1 We determine both the spin-modulation 
amplitude and decay length from the vortex core in fields up to H = 30 T. 



PACS numbers: 67.30.hm, 67.30.ht, 43.35.+d, 81.05.Rm 

Superconductivity and magnetism have been in the 
forefront of the study of high temperature supercon- 
ductivity, largely because their interplay is central to 
the pairing interaction from which superconductivity 
originates. An important manifestation of this di- 
chotomy can be found near a vortex core where cir- 
culating supercurrents create a sufficiently large mag- 
netic field as to suppress the superconducting state, 
thereby tipping the balance for stability in favor of 
magnetism. (2| A number of spin-sensitive, spatially- 
resolved, probes have been employed to explore this 
vortex state, including muon spin resonance, small an- 
gle neutron scattering, and nuclear magnetic resonance 
(NMR). These experiments can give the spatial distri- 
bution of local magnetic fields, as in the example of the 
17 NMR spectrum shown in Fig. 1. This spectrum, as 
well as results of other methods, are a superposition of 
local fields from both magnetism and superconductiv- 
ity and their deconvolution can be quite difficult. 

Here, we develop a new approach based on NMR 
spin-lattice relaxation (SLR) rate spectroscopy. The 
source for relaxation is the Doppler effect on the nodal 
quasiparticle excitations and is given by the square of 
the supercurrent momentum, p 2 s . Consequently, SLR 
can be attributed solely to superconductivity. Being 
two-dimensional, the rate is fully resolved from one 
superconducting plane to another, a decided benefit in 



comparison with measurements of the local field which 
average in three dimensions. Most importantly, SLR 
has strong spatial contrast increasing with the inverse 
square distance from the vortex core, ~ r~ 2 . Taking 
advantage of this spectroscopy, together with measure- 
ments of the NMR spectrum of local magnetic fields, 
we are able to characterize the magnetic spin den- 
sity wave (SDW) associated with vortices in our sam- 
ple, consisting of a slightly over-doped, single crys- 
tal of the anisotropic high temperature superconductor, 
Bi 2 Sr 2 CaCu20 8+ «5 (Bi2212), T c = 82 K. Our mea- 
surements extend over a wide range of applied mag- 
netic field, 4 T < H < 30 T. 

Scanning tunneling microscopy (STM) experiments 
on Bi2212 at moderate magnetic fields and low temper- 
atures showed that vortices induce a checkerboard pat- 
tern of the local density of states with a spatial period of 
15 A ss 4ao.|H| Further STM studies found such elec- 
tronic ordering in Bi2201[4j and NaCCOCE) persist- 
ing to non-superconducting dopings and above the su- 
perconducting transition temperature. Elastic neutron 
scattering experiments|6H8) on the cuprate supercon- 
ductor, La214, produced magnetic Bragg reflections 
corresponding to a periodicity of 30 Aw 8ao- Although 
spatially unresolved, this result was associated with 
the checkerboard pattern found from STM on Bi2212, 
since the period of a spin density wave is expected to 
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FIG. 1: The local field distribution and relaxation rate distribution in the vortex unit cell at H — 14.9 T, (a) calculated 
from Ginzburg-Landau (GL) theory. 1 3 1 The super current momentum squared is proportional to the Doppler term, and is the 
dominant contribution, in the spin-lattice relaxation rate, (b) The relaxation rates across the 17 NMR spectrum at T = 4.2 K 
of a slightly overdoped, T c = 82 K Bi2212 single crystal. On the low frequency side of the peak in the spectrum there is an 
unusual and substantial increase in the spin-lattice relaxation rate, T^ 1 , mirrored in the spin-spin relaxation rate, T 2 _1 . 



be twice that found in the density of states. 

There have been many theoretical interpretations 
for these results including competing SDW and super- 
conducting order[9 10], a Wigner crystal of Cooper 
pairs/pair density wave ATI [121 . striped ordering! 13 1, 
and d-density wave checkerboard order[14|. More re- 
cently, a phenomenological octet model of quasiparti- 
cle scattering from the fermi arcs [ 15 1 has gained pop- 
ularity in interpreting STM and STS results.! 16| The 
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FIG. 2: The NMR linewidth as a function of applied mag- 
netic field at 4.2 K. The contribution to the local field dis- 
tribution from supercurrents was calculated from GL theory 
(solid curve) and was subtracted from the data points. 



octet model fits some STM, STS, and ARPESlfTTl re- 
sults, however, this model does not account for the neu- 
tron scattering results[7| as well as the STM results 
of NaCCOC0. It has been suggestedHSl that, by 
increasing the applied magnetic field, magnetic order 
might be "tuned" and further suppress superconduct- 
ing order. Since NMR is a bulk, spin-sensitive probe 
and can be performed up to high field, it is ideal for 
investigation of the vortex states of high temperature 
superconductors. In the present work, we exploit these 
capabilities of NMR in addition to its potential for spa- 
tial resolution. 

In the case of an ideal, non-magnetic, supercon- 
ductor with straight vortex lines that form a two- 
dimensional lattice, the lowest frequencies of the NMR 
spectrum, Fig. l(a,b), come from the resonant nuclei in 
the sample positioned furthest from the vortex cores; 
the highest frequency components are from the nuclei 
closest to the vortex cores where the field is highest; 
and the peak corresponds to a saddle-point in the distri- 
bution of frequencies u = -fH, where 7 is the nuclear 
gyromagnetic ratio. This relation between spatial po- 
sitions in the vortex lattice unit cell and the NMR fre- 
quency is well-established. [ 19] However, contributions 
to the local fields from magnetism associated with vor- 
tex cores must be accounted for independently. 

A clear indication that such contributions exist in 
Bi2212 is evident from the spectrum of the central tran- 
sition for the O(l) planar site of 17 NMR, Fig. 1(b) 
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FIG. 3: (a) The spatial distribution of the local magnetic fields near a vortex within a C11O2 plane calculated for: (top) a spin 
density wave (SDW), Eq. 1, and (bottom) vortex supercurrents from GL theory|3] for which the amplitude Adia has been 
increased by a factor of 3 with respect to Asdw f° r clarity. (b)The amplitude Asdw (squares) and decay length r (circles) 
for the best fit spectrum generated from the local magnetic fields in (a). There is an increase in Asdw and r nearly constant 
with increasing magnetic field. [Inset] A comparison of the experimental (black), GL solution (blue, dotted), and SDW+GL 
(dashed-red) spectra centered at the central peaks, H = 15 T. 



and in the linewidths, Fig. 2. We find the distribu- 
tion of local magnetic fields to be much greater than 
that expected for vortex supercurrents alone. More- 
over, the linewidth increases systematically with ap- 
plied magnetic field, rather than decreasing as would 
have been the case for supercurrents represented by the 
solid curve in Fig. 2. Secondly, in Fig. 1(b) and Fig. 4, 
the SLR rate, T{" , is strongly non-monotonic across 
the spectrum contrary to predictions ll20l for d-wave 
superconductors without magnetic vortices, Fig. 1(a). 
Since SLR is from Doppler-shifted, nodal quasiparti- 
cles, it must increase smoothly and monotonically with 
decreasing distance approaching the vortex core (see 
supporting online materials). Thirdly, it is expected 
that the dominant interaction for spin-spin relaxation 
is from vortex vibrations which are also determined by 
vortex supercurrents. [21 1 Although there is not yet a 
quantitative theory for spin-spin relaxation, nonethe- 
less it should mimic SLR. From Fig. 1(b), the strong 
similarity between the two rate profiles indicates that 
this is correct. As with SLR, the spin-spin relaxation 
rate is not a monotonic function of local field. 

From these results we infer that the local field dis- 
tribution is not that of an ideal superconductor and has 
substantial magnetic contributions. Furthermore, the 
well-defined correlation between NMR relaxation and 
position in the spectrum requires that these additional 



contributions be closely associated with vortices. In 
fact, the non-monotonic behavior of the relaxation rate 
profiles can only be understood if the vortex local fields 
are oscillatory. We find that a spin-density wave model 
for vortex magnetism can completely describe both 
field distribution (spectrum) and supercurrent distribu- 
tion (SLR) and thereby characterize the vortex, spin- 
density wave state in Bi2212. 

In order to model the measured local field distribu- 
tions, we construct a spin "checkerboard" cosine wave 
decaying as a gaussian away from the vortex core. 
Consistent with STM and elastic neutron scattering 
measurements! 1 , 7| we take the period to be A = 30 
A~ 8ao and allow the amplitude and range of the spin- 
density modulation to be free parameters for fitting to 
the NMR spectra, 

B(x,y) = A SDW cos(2Trx/\)cos(2iry/\)e ( -- ( - x2+y2 V 2T2 '> 

(1) 

where Asdw ls the amplitude of the magnetization at 
the vortex core and r is its decay length. Eq. 1 pro- 
duces a SDW magnetic field component as shown in 
the upper part of Fig. 3(a). 

We also include in our fitting the local magnetic 
field distribution from diamagnetism of the vortex su- 
percurrents and we add this contribution to Eq. 1 at 
each position in the vortex unit cell before fitting to 
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FIG. 4: The calculated (red) and experimental (black) results for the spectra (curves) and relaxation rates (dots, model; squares, 
experiment) for various H at 4.2 K, (a-d) for (H — 4, 10, 15, 25 T). At low magnetic fields there is a near monotonic relation- 
ship for T^ 1 as a function of frequency (local field) as expected for the case Asdw < Ad ia , see Fig. 1(a). At high magnetic 
fields the SDW dominates producing a non-monotonic distribution of relaxation rates. It is clear that structure at high frequency 
in the experimental rate is replicated in the model calculation. This structure is sensitive to the spin wave period (see supporting 
online materials). The spectra are normalized to peak height and the model calculation for T^ 1 is offset vertically and scaled 
to compare with experiment. 



the data. The supercurrent diamagnetism was calcu- 
lated using Brandt's algorithm obtained by solving the 
Ginzburg-Landau (GL) equations for an ideal vortex 
lattice|3 1 and is shown at 15 T in Fig. 3(a) for a quanti- 
tative comparison with the corresponding SDW. To find 
the proper fitting parameters, we calculated the spectra 
for the GL solution plus SDW for many different values 
of Asdw an d T an d we used a chi-squared test to find 
the parameters shown in Fig. 3(b). With increasing ap- 
plied field H there is an increase in Asdw an d a nearly 
constant r, approximately twice the superconducting 
coherence length, r « 2£ Q sa 22 A. We then calculated 
the SLR rate being proportional to p 2 s , plus a constant 
offset, at each applied magnetic field where the super- 
current momentum is also calculated from GL theory. 
Multiple values of the calculated SLR rate, within each 
of 32 contiguous intervals of local field, were averaged 
and then presented as a single value to be compared 



with the data in Fig. 4. Over the entire range of applied 
field we capture an excellent representation of both the 
field distribution and the relaxation profiles with just 
two parameters, despite the unusual non-monotonic be- 
havior of the SLR rate as a function of frequency (local 
field). 

There are three ranges of applied magnetic field 
which characterize different profiles of T^ 1 across 
the NMR spectrum, depending on the relative contri- 
butions of the diamagnetic local fields Adia and the 
vortex-induced magnetic contribution Asdw ■ At low 
magnetic fields, 4-6 T, when Asdw < Adi a , re- 
sembles the near monotonic behavior expected for a 
e?-wave superconductor without magnetic vortices, 1 20 1 
see Fig. 1(a). Upon increasing the external magnetic 
field to 6-12 T there is a sharp increase of T± on the 
low frequency side of the spectrum in the range where 
Asdw ^ A dia . Finally, at high magnetic fields, 15-30 
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T, there is a significant increase of Tj~ on both low 
and high frequency sides of the spectrum correspond- 
ing to Asdw » Adia- Measurements in YBCO have 
revealed a small upturn in T^ 1 at low frequency, but 
much less pronounced than what we report here for 
Bi2212. In the YBCO work, contamination on the low 
frequency side of the spectrum from quadrupolar tran- 
sitions precluded an unobstructed view of the magnetic 
field distributions over the same wide range as is possi- 
ble with Bi2212. We think it is possible that vortices in 
YBCO are accompanied by a SDW checkerboard pat- 
tern which has not been resolved in these earlier NMR 
experiments. 

In summary we have shown that a spin-density wave 
can be associated with vortices in Bi2212 which ac- 
counts for both the NMR spectrum and NMR relax- 
ation behavior as a function of local magnetic field. 
The SDW amplitude increases with applied field and 
its spatial range is nearly constant at m 2£ Q up to 
H = 30 T. These results establish a link between 
previously reported checkerboard STM patternsfT| and 
spin-dependent elastic neutron scattering in a cuprate 
superconductor! 7] an d provide a detailed characteriza- 
tion of magnetic vortices in Bi2212. 
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SUPPORTING ONLINE MATERIALS 
MATERIALS AND METHODS 

Our Bi2Sr2CaCu20s+5 sample was post-processed 
by isotope exchange in 1 bar of flowing 17 (« 90% 
enriched) at 600 °C for 48 h followed by annealing for 
150 h at 450 °C resulting in a superconducting transi- 
tion temperature T c = 82 K. Spatially resolved T\ and 
T2 measurements were taken by Fourier transforming 
progressive saturation and Hahn echo sequences, re- 
spectively, with 16-pulse, phase alternation. Spectra 
from the central transition of the planar oxygen O(l) 
were acquired for a wide range of relaxation recov- 
ery times and then the frequency domain of each one 
was subdivided into 32 intervals. The recovery pro- 
files in each interval were fit to a progressive saturation 
recovery J22] for Ti and an exponential for T 2 . Exper- 
iments presented here were performed at the National 
High Magnetic Field Lab in Tallahassee, Florida and 
at Northwestern University for fields from 4 to 30 T at 
a temperature of 4.2 K. The spectra were taken using 
a frequency sweep technique. Before taking data at a 
particular applied magnetic field, we heated the sam- 
ple above T c for 10 min followed by quench cooling 



to 4.2 K. This procedure consistently produced repro- 
ducibly narrow spectra through many field cycles for 
which one indication is the systematic dependence of 
the linewidth on magnetic field, Fig. 2. Data for some 
fields were also acquired at different temperatures, but 
are not reported here. Measurements were also taken 
at magnetic fields below those shown in Fig. 2 and are 
reported separately 11231 . 

Spectra from the oxygen sites, [24] O(l) in the CuC>2 
plane and 0(2) in the SrO plane, are well-resolved at 
T > 40 K. At lower temperatures the central transi- 
tions overlap. Discrimination between O(l) and 0(2) 
can be achieved by saturation of the 0(2) resonance 
which is only weakly coupled to the electronic excita- 
tions in the CuC>2 plane and has an order of magnitude 
smaller compared to O(l). The relative intensi- 
ties of the central and quadrupolar satellites in the O(l) 
spectrum, which do not overlap with the 0(2) spec- 
trum, are monitored down to 4.2 K to assure that the 
O(l) central transition has no contribution from 0(2). 

NMR LINEWIDTH 

The temperature of the transition from vortex liquid 
to vortex solid has been determined as a function of 
magnetic field by Chen et al. [24| in Bi2212, signaled 
by an abrupt increase in linewidth with decreasing tem- 
perature. In the liquid state the time scale for vortex dy- 
namics is much shorter than that of the NMR measure- 
ment, thereby averaging the local magnetic fields asso- 
ciated with vortices. The experiments reported here are 
performed in the vortex solid phase. 

It is instructive to compare the linewidth for 17 
NMR reported here with that measured in YBCO in the 
normal state for various cations (Zn, Ni, and Li) substi- 
tuted for copper in the CuC<2 plane generating an impu- 
rity local moment that contributes to the lineshape. The 
17 O NMR broadening we observe in Bi2212 is similar 
in magnitude to that reported for Zn in YBCO at the 
same field. For this comparison we make an adjust- 
ment for the vortex concentration (Bi2212) to be the 
same as the cation impurity concentration (YBCO). We 
assume the same hyperfine coupling for Bi2212 and in- 
fer, based on the NMR spectrum alone, that the vortex 
core could be modeled as a local impurity with moment 
of ~ 0.6/is. By local impurity we mean a perturbation 
that induces local field, RKKY-type oscillations with 
period ~ ir/kp, which is about a factor of four smaller 
than the modulation period we use here. Although our 
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analysis does not allow the numerical freedom to in- 
dependently choose the period of the modulation for 
all applied fields, nonetheless, at the highest magnetic 
fields we can perform a x 2 analysis to simultaneously 
fit the measured SLR rate and the spectrum in terms of 
the model of Eq. 1. We find that the oscillation period 
is then constrained to be 30 ± 5 A, as might be ex- 
pected from the reports from STM and elastic neutron 
scattering. (U |7] 

NMR RELAXATION 

Spatially resolved NMR relaxation experiments 
have been carried out on YBCO aligned powders in 
previous work. [22, 26, 27 1 However, the present single 
crystal samples of BSCCO show significantly higher 
quality results, in part because the quadrupolar satel- 
lites are fully resolved. Nonetheless, both the YBCO 
and Bi2212 measurements at low temperatures, T < 
15 K, and at high magnetic fields, 10 < H < 30 
T, consistently indicate ll22l l27l that at each frequency 
position in the spectrum, T\T is a constant and that 
Ti increases with applied magnetic field uniformly 
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across the spectrum, except in the vortex core re- 
gion. This quadratic field dependence has been at- 
tributed to the Zeeman shift in the nodal quasiparti- 
cle (NQ) energy. Il22ll27l The mechanism for relaxation 
with nodal quasiparticles includes thermal, Zeeman, 
and Doppler contributions. Based on the above obser- 
vations, at fixed applied magnetic field, the dependence 
of the SLR rate on position in the vortex unit cell can be 
attributed to the Doppler shift in the NQ energy, vj?.p s . 
Since the SLR rate depends on the square of the local 
density of states (LDOS) and the LDOS is linear in the 
NQ energy, then the Doppler contribution to the SLR 
rate is cx p 2 as displayed in Fig. 1(a). These arguments 
are presented in more detail by Mitrovic et al.[22\ 

The time derivative of the measured relaxation pro- 
file in the limit of short times, averaged over the entire 
spectrum, gives the average SLR rate. This is propor- 
tional to the average of the square of the LDOS over the 
vortex unit cell whose spatial area is itself proportional 
to applied magnetic field. Volovik[25| pointed out that 
such an average is dominated by the Doppler term in 
the electronic excitation energy. In Fig. 5 we show the 
magnetic field dependence of the average SLR rate, red 
points, together with the Zeeman contribution to the 
rate, green points, measured at the peak in the spec- 
trum, i.e. the saddle point. The difference between the 
data sets is due to the Volovik effect following the ex- 
pected H logiJ behavior which we show, together with 
the Zeeman contribution, as the (red) dashed curve in 
the figure. This agreement confirms that the Doppler 
effect on NQ is a principal mechanism for relaxation 
and depends on position in the vortex unit cell. Con- 
sequently, a measurement of can distinguish be- 
tween atoms in the NMR signal that are at different 
positions relative to the vortex core thereby providing 
spatial resolution within the CUO2 plane. 



FIG. S 1 : The average SLR rate, T{~ , as a function of mag- 
netic field. The average over the vortex unit cell (red points) 
was obtained by weighting the measured T{~ by the spec- 
trum, i.e. the probability distribution of local fields. We 
expect that the contribution to the average T{~ from the 
Volovik effect !25l is proportional to HlogH which, when 
added to the Zeeman contribution (green points, measured 
at the saddle point peak in the spectrum), fits the total av- 
erage field dependence as shown by the dashed curve. The 
intercept at H = is the thermal contribution to T^ 1 . 



